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Abstract We revisit Voloshin's model of multiple black hole production in trans- 
Planckian elementary particle collisions in D = 4. Our revised compu- 
tation shows that the cross section to produce n additional black holes is 
suppressed by s _1 , rather than being enhanced as was originally found. 
We also review the semiclassical gravity picture of black hole production 
from hep-th/0409131 making additional comments about the meaning 
of wavepacket subdivision. 



1. INTRODUCTION 

Black hole (BH) production in trans-Planckian elementary particle 
collisions (E S> Epi anc ^) has long been considered a theoretical possi- 
bility. If TeV-scale gravity scenarios based on large extra dimensions or 
warped compactifications are realized in nature, this possibility may be 
realized in practice at future accelerators (see tl for a recent review). 
The key question is the cross section of this process, which is usually 
assumed to be set by the horizon radius of the produced BH (the so 
called "geometric cross section" a ~ nrfy. 

In this note I would like to discuss two aspects of the BH produc- 
tion problem relevant for justifying the geometric cross section estimate. 
First I will explain how to derive this estimate from a controlled semi- 
classical gravity approximation to the BH production process, adding 
some comments to the original discussion of [2] (also recently reviewed 
in [3]). Then I will estimate the cross section of multiple BH production 
due to collisions between virtual gravitons emitted by the primary parti- 
cles, in a model first proposed by Voloshin & i . The conclusion (contrary 
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to is that multiple BH production gives a subdominant contribution 
to the total cross section. 

The size of produced BHs in large extra dimension scenarios is typi- 
cally much smaller than the size of extra dimensions, and their produc- 
tion may thus be considered as happening in flat D dimensional space- 
time. Since we are focussing on theoretical issues, to keep the discussion 
clear we will work in D = 4. We will also use Planck units, setting 

E ' Planck = 1- 

2. WAVEPACKETS AND SEMICLASSICS 

In this note we will adhere to the standard believe that in quantum 
theory of gravity classical BHs with mass M S> 1 will be realized as long- 
lived resonance states, decaying via Hawking radiation, with lifetime 
~ M 3 . It is then energetically allowed to produce such BHs in trans- 
Planckian {E 1) elementary particle collisions. 

To estimate the cross section of this process, Eardley and Giddings 5 
looked at the grazing collision of two ultrarelativistic point particles in 
classical general relativity, using formation of a closed trapped surface 
(CTS) as a sufficient condition for BH formation. In this totally classical 
description, lower bound for the cross section is given by ^h 2 max ^ where 
bmax is the maximal impact parameter for which we are able to find 
a CTS in the spacetime formed by two colliding Aichelburg-Sexl shock 
waves. [3] found b max ~ E ~ r^, and this implies a ~ itr\. Recourse 
to such an indirect method is necessary, because explicit solutions of 
Einstein's equations exhibiting the final BH state are out of reach. 

How do we justify this approach from quantum gravity point of view? 
A point of immediate concern is that particles in a collider experiment 
are described by wide wavepacket states of macroscopic size (set essen- 
tially by the beam radius, which in turn is determined by the focussing 
ability of accelerator magnets). These wavepackets are vastly larger 
than BHs whose production we are trying to describe. If we use the 
energy momentum tensor of these wavepacket states \ip) in the RHS of 
the semiclassical Einstein's equations 

V- lg fi vR = &\T^ (i.i) 

we won't see any BH production whatsoever, since the energy is spread 
out over a huge volume, and the energy density is insufficient to cause 
collapse. Informally speaking, particles "do not fit" inside a BH. 

However, this does not mean that BHs do not form. The correct in- 
terpretation is that the part of the gravitational field wavefunction cor- 
responding to BH production was erased — averaged away — by eq. (jl.lj) . 
Some averaging is always inevitable when using the semiclassical field 
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Figure 1.1 Subdivision of the right-moving particle wavepacket. Similar subdivision 
has to be done for the left-moving particle. 



equations, since we substitute T^ v by its expectation value. Unfortu- 
nately, in this case it destroys precisely the part of the wavefunction we 
are interested in. 

To see BH production, one should instead proceed as follows. First 
of all, we have to subdivide the initial wavepackets into much smaller 
wavepackets of size w <C r^: 

|^)=iV- 1 / 2 ^ |^>. (1.2) 

This subdivision is carried out so that the small wavepackets \ipi) in 
the RHS are almost orthogonal. This orthogonality is quite obvious in 
the position representation (see Fig. II. lj) . Because of the orthogonal- 
ity, collisions between different pairs of small wavepackets are mutually 
excluded possibilities, and probabilities of BH production in each such 
elementary collision should be added. Now, it is the collisions of the 
small wavepackets that we are going to analyze using eq. (jl.ljl . Con- 
dition w <C r/j ensures that the small wavepackets produce a collision 
spacetime which is a small perturbation of the one corresponding to 
point particles of the same energy. Thus the Eardley-Giddings analysis 
applies, and adding probabilities results in the geometric cross section. 

Using finite-size wavepackets instead of point particles has an addi- 
tional bonus in that it puts the conditions of applicability of the semi- 
classical approximation under control. For example, curvature blows up 
when the shock fronts of Aichelburg-Sexl waves corresponding to point 
particles collide [HI Ej. However, taking wavepacket size into account 
regulates the curvature and brings it below the Planck value, so that we 
can trust the semiclassical gravity approximation |2j. 

3. MULTIPLE BLACK HOLE PRODUCTION 

According to the above discussion, the geometric cross section formula 
provides a lower bound for a single large BH production cross section 
in a trans-Planckian collision. However, as the energy of the particles 
grows, multiple BH production also becomes energetically allowed. It is 
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Figure 1.2 A typical diagram for multiple BH production in Voloshin's model 

important to understand which process is dominant at asymptotically 
high energies. If multiple BHs dominate, it will be much harder to 
observe dipole patterns of emitted particles expected in the Hawking 
evaporation of a single large BH. 

We are going to discuss a model proposed by Voloshin |3j, in which 
multiple BHs are produced due to collisions between virtual gravitons 
emitted by the trans-Plan ckian projectiles. Such virtual graviton emis- 
sion is a quantum effect: classically, any radiation happens after the 
particles collide. The process is studied diagrammatically, with a typ- 
ical diagram shown in Fig. 11.21 Only the case of small peripheral BHs 
(1 <C rrii <C E, i = 1 . . . n) is considered, so that the gravitons are "soft". 

The production amplitude is computed from the diagrams using the 
standard QFT propagators and particular vertices for soft graviton emis- 
sion and for BH production (Fig. I1.3jl . The amplitude to emit a posi- 
tive helicity graviton of energy uj <C E and small transverse momentum 
k = (/c2, ^3), |k| <C u> is given by (see Appendix A) 



(We will not pay attention to constant numerical factors. Thus our final 
result Q1.5|) is valid up to a factor c n .) 

Although we are unable to compute the elementary BH production 
vertex f(q 2 ), the geometric cross section allows to fix the combination 



A crucial final element of the model is a condition which ensures that 
the emitted gravitons do not subsequently fall into a common large BH. 



Aoc {E/u>) 2 (k 2 + ik 3 ) 2 . 



(1.3) 




(1.4) 




Figure 1.3 Graviton emission and BH production vertices 
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Such an infall may happen due to graviton rescattering diagrams, which 
we are not going to compute. Thus, without a "fall safe" condition we 
would be in danger of greatly overestimating the multiple BH amplitude. 
Voloshin's "fall safe" condition limits the transverse momenta of emitted 
gravitons by |k| < 1/E. 

To derive this condition, we note that a typical emitted graviton will 
be off-shell by AE ~ k 2 /u;. It will exist for a time interval At ~ 1/AE, 
during which it will reach transverse separation Az ~ (k/uj)At from the 
projectile. Voloshin's condition arises if we require that this transverse 
separation is larger than the horizon radius of the big BH formed in the 
collision of the primary particles: Az > E. 

In the described model, our computation (see Appendix B) gives the 
following amplitude to produce n additional BHs with 4- momenta qf. 

f in \s,qf) ~ f(s)H n i=i (qf)- 2 f(qf) . (1.5) 

The original computation of jlj, Eq. (5), gave an amplitude larger than 
(|1.5|) by a factor of (sqf) n . We believe that our result is correct; see 
Appendix B and [£>] for an explanation. 

Using Q1.4JI . we can compute from ()1.5j) the contribution of the dia- 
grams from Fig. 11.21 into the total cross section a n to produce one large 
(m 2 ~ s) and n small BHs. This contribution will behave like s l ~ n , the 
suppression being due to the phase space restriction |k| < 1/E satisfied 
by the small BHs as a consequence of the "fall safe" condition. 

However, for n > 2 there are diagrams which give a larger contri- 
bution, so that <7 n ~ const is likely for any n > 1. Consider, e.g., 
Fig. 11.41 where the primary particles emit "fall safe" gravitons of en- 
ergy Ei, 1 <C E\ <C E, and it is these gravitons that form n smaller 
BHs according to the previous model. The allowed phase space for this 
diagram will be much bigger, since the individual small BHs can now 
have much larger transverse momenta |k| < 1/E\, only their sum being 
< 1/E. Choosing Ei above the threshold of n BH production, we will 
get an s-independent contribution to an. 1 

In any case, we see that a n is suppressed compared to the geometric 
cross section value gq ~ s. 



1 It is also easy to see that the inclusive cross section of multiple BH production cannot decay 
with s. This is because the particles may first reduce their energy by emitting one or more 
"fall safe" gravitons (which costs no s-dependent factor, see [3], Eq. (3)), and then collide to 
form BHs. I am grateful to M. Voloshin for this remark. 
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Figure 1.4 A diagram giving an s-independent contribution to a n 



4. CONCLUSIONS 

In this note we discussed two aspects of the BH production problem. 
In Section 2 we discussed how wavepacket arguments can be used to 
justify the use of semiclassical gravity in this problem. In Section 3, we 
revisited Voloshin's model of multiple BH production, and presented a 
revised computation which shows that this process is suppressed com- 
pared to the production of a single large BH. Our main conclusion is 
that the geometric cross section estimate seems to be in rather good 
health, surviving all checks and resisting any disproving attempts. 
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Appendix: A. Soft graviton emission vertex 

Consider emission of a positive helicity collinear graviton with energy 
lj and small transverse momentum k = (&2,0), k 2 <C uj. Its polarization 
tensor, satisfying the constraints hijkj = 0, ha = 0, is given by 

(k^/uj 2 —kijio —ik2/u) \ 
-hi/ui 1 i (l.A.l) 

-ik 2 /uJ i -\-k\ju) 2 ) 

The emission amplitude is A oc (p — k\Tij\p)hij, where the energy mo- 
mentum tensor matrix element is 

(P\T^\p) = (P|0>, v -i^(90) 2 |p) =PvJ> v +PvPp-V^(p p )- ( LA - 2 ) 
From this we find A oc E 2 k\juP'. 
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Appendix: B. Multiple BH production amplitude 

We will compute the n = 2 amplitude in the case when both small BHs 
are produced at rest in the cm. frame: qi = (mj, 0,0,0) in Fig. II .B. II 
The calculation of the amplitude naturally splits into 3 steps. First, 




Figure l.B.l A diagram for one large and two small BHs in the final state 



we compute the loop integral over the longitudinal momenta. Then, 
we multiply by the emission vertices and integrate over the transverse 
momenta. Finally, we multiply by the BH production vertices. 

It is convenient to write the graviton momenta as follows {i = 1, 2): 

ki = (^ L + ^Y t ,^--^Y t M, K = qi -h. (i.b.i) 

The longitudinal loop integral now separates into the integrals over xi 
and yi. The part depending on Xi is 

/= /^ 2 ^(p-w4>-*i-fe) 2 (LB ' 2) 



dx\ dx2 



(misi - kf)(m 2 x 2 - kl)(-Ex 1 - kf)(-E(xi + x 2 ) - kf 2 ) 



The +ie is implicit in each denominator. This integral is easy to compute 
by closing the contour in the lower half-plane. Omitting corrections of 
the order m/E <^ 1, we have 

I cx^kjikj + kl)]- 1 . (1.B.3) 

Since we used complex analysis, it is important to check that the 
integral is dominated by soft, almost real gravitons. This is indeed true, 
the important region being \x{\ < k?/mj. The tail corresponding to 
\x.- L \ S> kf/mj can be estimated directly as: 

1 f dxi dx2 T „ .v 

' -o (1.B.4) 
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This check also justifies post factum neglecting the dependence of the 
graviton emission vertices on Xj, as well as the omission of x 2 terms in 
the denominators of (jl.B.2|) . 

Before proceeding to the next step, we have to add diagrams differing 
by the order of graviton emission (Fig. I1.B~2|) . This summation has the 




Figure 1.B.2 One of the 3 permuted diagrams 



effect I — > [-E 2 k^k|] . Finally, we multiply by the same factor arising 
from the y- integration, and get: 

[E 2 \4\4]- 2 (longitudinal part) . (1.B.5) 

We will assume that the small BHs are produced in the spin state, 
so that the colliding gravitons have opposite helicities. The product of 
the corresponding emission vertices is 

oc {E/m)\k 2 + ik 3 ) 2 (k 2 - ik 3 ) 2 = (E/m)\k 2 ) 2 . (1.B.6) 

Multiplying H1.B,5|) by such factors for both BHs and integrating over 
[kj| < 1/E (the "fall safe" condition), we get a number ~ (mim2)" 4 . 
Finally, we multiply by the BH production vertices and arrive at 

/ (2) ~ f(s) f(m 2 ) f(m 2 ) (m im2 )- 4 . (1.B.7) 

This formula agrees with the general result (|1.5|) in the considered case. 

Extension of the above computation to the general case is quite straight- 
forward. A point worth mentioning is the use of the standard identity 

V [ ai ( ai + a 2 ) ■ ■ ■ {a! + . . . + an)}" 1 = [ ai ■ ■ ■ an}' 1 (1.B.8) 

<—^perm 

when summing over the order of graviton emission in the n > 2 case. 

In m, the estimate J d ki/(kfk' 2 ) ~ 0(1) was used in computing the 
amplitude. However, our analysis (see also |S]) shows that the correct 
estimate is: 

d 4 k f dxdyd 2 k 1 

<* / 7ZIZ 1.. I 7±7—. I 73 ~ • (1-B.9 



k 2 k' 2 J (mx — k 2 + ie)(my — k 2 + ie) m 2 E 2 

This extra factor, for each of n BHs, explains the difference between 
lfT5Jl and the result of gj. 
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